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ABSTRACT 

We present adaptive optics assisted integral field spectroscopy of nine Ha-selected galaxies at 
z = 0.84-2.23 drawn from the HiZELS narrow-band survey. Our observations map the kinematics 
of these star-forming galaxies on ~ kpc-scales. We demonstrate that within the ISM of these galaxies, 

the velocity dispersion of the star-forming gas (a) follows a scaling relation a oc Sgpp + constant 
(where Esfr is the star formation surface density and the constant includes the stellar sur- 
face density). Assuming the disks are marginally stable (Toomre Q — l), this follows from the 
Kennicutt-Schmidt relation (SgpR = ASg l as ), and we derive best fit parameters of n— 1.34±0.15 and 

A = 3.4ti'e x 10~ 4 Mq yr _1 kpc~ 2 , consistent with the local relation, and implying cold molecular gas 
masses of M gas = 10 9 ~ 10 M Q and molecular gas fractions M gas / (M gas + M*) = 0.3 ±0.1, with a range 
of 10-75%. We also identify eleven ~kpc-scale star-forming regions (clumps) within our sample and 
show that their sizes are comparable to the wavelength of the fastest growing mode. The luminosities 
and velocity dispersions of these clumps follow the same scaling relations as local Hn regions, although 
their star formation densities are a factor ~ 15 ± 5 x higher than typically found locally. We discuss 
how the clump properties are related to the disk, and show that their high masses and luminosities 
are a consequence of the high disk surface density. 

Subject headings: galaxies: starburst, galaxies: evolution, galaxies: high-redshift 



1. INTRODUCTION 



The majority of 
sive galaxies (M* 
years ago, 



the stars 
>1O U M ) 
an epoch 



in the most mas- 
formed around 8 - 
10 billion years ago, an epoch when star formation 
was at its peak (|Hopkins fc Beaco~ml 120061 : iSobral et al.l 
l2012bf) . Galaxies at this epoch appea r to be gas- 



rich ( L as = 20 - 80%: iTacconi et al.l l20lcl iDaddi et al 



20101: [Geach et all 120 111) and turbulent ([Lehnert eta! 



2009(1 . with high velocity dispersions given their ro- 
tation al velocities (a = 30- 100 km s" 1 , iWx la ~0-2- 
1: e.g. iForster Schreiber et al.l 120091: IGenzel et al.l I2008t 
IWisnioski et alj|2011t IBothwell et al.ll2012D . Within the 
dense and highly pressurised inter-stellar medium (ISM) 
of these high-redshift galaxies, it has been suggestes that 
star formation may be triggered by fragmentation of dy- 
namically unstable gas (in contrast to star-formation oc- 
curring in giant molecular clouds in the Milky- Way which 
continually condense from a stable disk and then dis- 
sipate). This process may lead to the to the forma- 
tion of massive (~10 8 ~ 9 M(T)) star-forming regions (e.g . 
lElmegreen et al.l 120071: IBournaud fe Elmegreenl 120091) 
and give rise to the the clumpy morp hologies that are 
often seen in high-redshift starbursts (Elmcgrecn et al.l 
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l200l . 

In order to explain the ubiquity of "clumpy" disks 
seen in images of high-redshift galaxies, numerical sim- 
ulations have also suggested that most massive, star- 
forming galaxies at z = 1- 3 continually accrete gas from 
the inter-galactic me dium along cold and clumpy streams 
from the cosmic web (IKeres et al.l2005HDekel et al 112 009: 
IBournaud fe Elmegreenll2009Uvan de Voort et all 2011). 
This mode of accretion is at its most efficient at z ~ 1 - 2, 
and offers a natural route for maintaining the high gas 
surface densities, star formation rates and clumpy mor- 
phologies of galaxies at these epochs. In such models, 
the gas disks fragment into a few bound clumps which 
are a factor 10- 100 x more massive than star-forming 
complexes in local galaxies. The gravitational release of 
energy as the most massive clumps form, torques between 
in-spiraling clumps and energy injection from star forma- 
tion are all likely to contribute to maintaining the high 
turbulence yelocity dispersion of the inter-stellar medium 
(IgM ) (e.g. IBournaud fc Ermegree n]|2009l: iLeh nert et al.l 
12001 iGehzel et al.ll2008L I2011D . 

In order to refine or refute these models, the obser- 
vational challenge is now to quantitatively measure the 
internal properties of high-redshift galaxies, such as their 
cold molecular gas mass and surface density, disk scal- 
ing relations, chemical make up, and distribution and 
intensity of star formation. Indeed, constraining the evo- 
lution of the star formation and gas scaling relations 
with redshift, stellar mass and/or gas fraction are re- 
quired in order to understand star formation throughout 
the Universe. In particular, such observations are vi- 
tal to determine if the prescriptions for star formation 
which have been developed at z = can be applied to 
the rapidly evolving ISM of gas-rich, high-redshift galax- 
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ies (|Krumholz fe Dekell[20Tot IHopkinsll20ll) . 

To gain a census of the dominant route by which galax- 
ies assemble the bulk of their stellar mass within a well 
selected sample of high-redshift galaxies, we have con- 
ducted a wide field (several degree-scale) near-infrared 
narrow-band survey (the High-Z Emission Line Survey; 
HiZELS) which targets Ha emitting galaxies in four pre- 
cise ( Az = 0.03) redshift slices: z = 0.40, 84, 1.47 and 
2.23 (IGeach et all l2008t ISobral et al.ll200i [20lol I20T1 
[2012a|bJ). This" survey provides a large, star formation 
limited sample of identically selected Ha emitters with 
properties "typical" of galaxies which will likely evolve 
into ~ L* galaxies by z = 0, but seen at a time when they 
are assembling the bulk of their stellar mass, and thus 
at a critical stage in their evolutionary history. More- 
over, since HiZELS was carried out in the best-studied 
extra-galactic survey fields, there is a wealth of multi- 
wavelength data, including 16-36 medium and broad- 
band photometry (from rest-frame UV - mid- infrared 
wavelengths allowing robust stellar masses to be de- 
rived), Herschel 250-500/im imaging (allowing bolomet- 
ric luminosities and star formation rates to be derived) 
as well as high-resolution morphologies for a subset from 
the Hubble Space Telescope CANDELS and COSMOS 
ACS surveys. 

In this paper, we present adaptive optics assisted inte- 
gral field spectroscopy of nine star-forming galaxies se- 
lected from HiZELS. The galaxies studied here have Ha- 
derived star formation rates of l-27M Q yr _1 and will 
likely evolve into ~ L* galaxies by z = 0. They are there- 
fore representative of the high-redshift star-forming pop- 
ulation. We use the data to explore the scaling relations 
between the star formation distribution intensity and gas 
dynamics within the ISM, as well as the properties of 
the largest star-forming regions. We adopt a cosmology 
with f2 A =0.73, fl m = 0.27, and H = 72 kms -1 Mpc -1 in 
which 0.12" corresponds to a physical scale of 0.8 kpc at 
z = 1.47, the median redshift of our survey. All quoted 
magnitudes are on the AB system. For all of the star 
formation rates and stellar mass estimates, we use a 
IChabrierl (|2003D initial mass function (IMF). 

2. OBSERVATIONS 

Details of the target selection, observ a tions and data- 
reduction are given in ISwinbank etal! ((2012). Briefly, 
we selected nine galaxies from HiZELS with Ha fluxes 
0.7- 1.6 x 10~ 16 ergs -1 cm~ 2 (star formation rate43of 
SFRfj Q = 1 — 27 Mq yr -1 ) which lie within 30" of bright 
(R<15) stars. We performed natural guide star adap- 
tive optics (AO) observations with the SINFONI IFU 
between 2009 September and 2011 April in ~ 0.6" seeing 
and photometric conditions with exposure times between 
3.6 to 13.4 ks. At the three redshift slices of our targets, 
z = 0.84[2], z = 1.47[6] and z = 2.23[l], the Ha emission 
line is redshifted to ~ 1.21, 1.61 and 2.12/xm (i.e. into 
the J, H and if-bands respectively). The median strehl 
achieved for our observations is 20% and the median en- 
circled energy within 0.1" (the approximate spatial res- 
olution of our observations) is 25%. 

The data were reduced using the SINFONI esorex 
data reduction pipeline which extracts, flat-fields, wave- 
length calibrates and forms the data-cube for each ex- 
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posure. The final (stacked) data-cube for each galaxy 
was generated by aligning the individual data-cubes and 
then combining them using an average with a 3-cr clip to 
reject cosmic rays. For flux calibration, standard stars 
were observed each night either immediately before or 
after the science exposures and were reduced in an iden- 
tical manner to the science observations. 

As Fig. [1] shows, all nine galaxies in our SINFONI- 
HiZELS survey (SHiZELS) display strong Ha emission, 
with luminosities of Ln a ~ 10 414_42 ' 4 ergs~ 1 . Fitting 
the Ha and [Nil] AA6548, 6583 emission lines pixel- by- 
pixel using a \ 2 minimisation procedure we construct 
intensity, velocity and velocity dispersi on maps of our 
samp le and show these in Fig. Q] (see also ISwinbank et all 
2012 for details). 

3. ANALYSIS & DISCUSSION 
3.1. Galaxy Dynamics and Star Formation 

As ISwinban k et al.1 ()2012l ) demonstrate, the ratio 
of dynamical-to-dispersion support for this sample is 
usin(i) / (7 = 0.3-3, with a median of 1.1±0.3, which is 
consistent with similar measurements for both AO and 
non- AO studies of star-forming ga laxies at this epoch 
(e.g. iForster Schreiber et all 120091 ) . The velocity fields 
and low kinemetry values of the SHiZELS galaxies (total 
velocity asymmetry, K to t =0.2-0.5) also suggest that at 
least six galaxies (SHiZELS 1, 7, 8, 9, 10, & 11) have 
dynamics consistent with large, rotating disks, although 
all display small-scale deviations from the best-fit dy- 
namical model, with <data— modcl> = 30±10kms _1 , 
with a range from <data — model> = 15-70kms _1 
(jSwinbank et al.ll2012t ). 

We also use the multi-wavelength imaging to 
calculate the rest-frame SEDs of the galaxies in 
our sample and so derive the stellar mass, red- 
dennin g and estimates of the star-formation his- 
tory (So bral et al.l 1201 If). From the broad-band 
SEDs (Fig. 1 of ISwinbank etUI I20T1 . the aver- 
age E(B-V) for our sample is E(B - V) = 0.28 ± 0.10 
which corresponds to A„ = l.ll±0.27mag and indi- 
cates An a = 0.91 ± 0.21 mag. The resulting dust- 
corrected Ha star formation rate for the sample is 
SFRh q = 16± 5Mq yr _1 , which is consistent with that 
inferred from the far-infrared SEDs using stacked Her- 
schel SPIRE obser vational SFR FT r = 18 ± 8 M yr -1 ; 
ISwinbank et al.ll2012j ) 

Next, to investigate the star formation occurring 
within the ISM of each galaxy, we measure the star for- 
mation surface density and velocity dispersion of each 
pixel in the maps. Since we do not have spatially re- 
solved reddening maps, for each galaxy we simply cor- 
rect the star formation rate in each pixel using the best- 
fit E(B — V) for that system. We also remove the ro- 
tational contribution to the line width at each pixel by 
calculating the local AV / AR across the point sp read 
function (PSF) for each pixel (jDavies eTffll2l)Tl . In 
Fig. [2] we plot the resulting line of sight velocity disper- 
sion (a) as a function of star formation surface density 
(Ssfr) for each galaxy in our sample. We see that there 
ap pears to be a corre l ation between Esfr and a, and 
as iKrumholz fe De kcl (2010) show, this power-law cor- 

4 Adopting SFRfir, = 2.7 x 1(T 44 L FIR (ergs- 1 ) 
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Fig. 1. — Ha intensity, velocity field, line of sight velocity dispersion (cr) and Toomre (Q) maps of the nine SHiZELS galaxies in our 
sample. Top Left: Ha emission line map. In SHiZELS 7, 8, 9, and 14 we identify and label the star-forming regions (clumps). Top Right: 
Toomre Q(x,y) maps of each galaxy, with contours at Q = 0.5, 1.0, and 2.0. In galaxies where we have identified star-forming regions 
(clumps), we also overlay their positions. These star-forming regions have an average Toomre Q = 0.8 ± 0.4. Bottom Left: Ha velocity field 
of each galaxy (with the best-fit kinematic model overlaid as contours). Bottom Right: line of sight velocity dispersion (ex), corrected for 
local velocity gradient (A V/ AR) across the PSF. At least six galaxies (SHiZELS 1, 7, 8, 9, 10, & 11), have dynamics that indicate that 
the ionised gas is in a large, rotating disk. A further two are compact (SHiZELS 4 & 12) whilst the dynamics of SHiZELS 14 indicate a 
merger. 



relation may be a natural consequence of the gas and 
star formation surface density scaling laws. For example, 
first consider the Toomre stability criterion, Q, (jToomrd 
[lQfil . 



Q = 



(T K 



TrGSd 



isk 



(i) 



where a denotes the line of sight velocity dispersion, Sdisk 
is the average surface density of the disk, k = a v max / R 
where w max is the rotational velocity of the disk, R is the 
disk radius and a = \pl for a flat rotation curve. Galaxy 
whose disks have Q < 1 are unstable to local gravita- 
tional collapse and will fragment into clumps, whereas 
those with Q Z 1 have sufficient rotatio nal support for 
the gas to wi thstand and collapse. As iHopkind (|2012j ) 
(e.g. see also iCacciato et al.l I2012D point out, gas-rich 



galaxies are usually driven to Q ~ 1 since regions with 
Q < 1 begin forming stars, leading to super-linear feed- 
back which eventually arrests further collapse due to 
energy/momentum injection (recovering Q ~ 1). For 
galaxies with Q ^>1, there is no collapse, no dense re- 
gions form and hence no star formation (and so such 
galaxies wo uld not b e selec ted as star-forming systems) . 

Following iRafikovl (|2001f ). and focusing on the largest 
unstable fluctuations, the appropriate combination of gas 
and stellar surface density (S gas and respectively) is 



J disk 



i + fs 



(2) 



where f a = <J*l<J g is the ratio of the velocity dispersion 
of the stellar component to that of th e gas (see also the 
discussion in lRomeo h WiegertJl2011|) .. 
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Next, iKennicuttl (|1998bT ) show that the gas and star 
formation surface densities follow a scaling relation 



J SFR 



M yr 1 kpc~ 



= A 



Jgas 
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Fig. 2. — Star formation rate surface density as a function of velocity dispersion for each pixel within the galaxies in our sample. The 
star format ion rates are derived from Ho, corrected for galaxy reddening and the velocity dispersion has been corrected for local velocity 
gradient (§ 13. ip . The small solid and open symbols denote measurements within and outside the half-light radius respectively. The solid 
squares show the star-formation and velocity dispersions of the ~ kpc-scale clumps (Table 2) which appear as regions of high star formation 
density given their velocity dispersion. The grey region denotes the best fit to the Esfr - t relation from combining the Toomre criterion 
and Kennicutt-Schmidt law (see equation (4} with power-law index ranging from n = 1.0- 1.4 (the dashed curve shows the solution for 
71 = 1.2). Over this range, the data is consistent with an absolute star formation efficiency of A = 4.1±2.4 X 10 — 4 Mg yr - 1 kpc -2 . 

In order to test whether this model provides an ad- 
equate description of our data, we fit the Esfr _ c 
distribution for each galaxy in our sample. To es- 
timate the s tellar surface density, E*, we we follow 
iSobral et all (poTTI ) and perform a full SED x 2 St of 
the rest-frame UV — mid- i nfrared photometry using the 
iBruzual fc Charlotl (f2003h and Bruzual (2007) popula- 
tion synthesis models. We use photometry from up to 
36 (COSMOS) and 16 (UDS) wide, medium and narrow 
bands (spanning GALEX far-UV and near-UV bands to 
Spitzer/IRAC) and calculate the rest-frame spectral en- 
ergy distributi on, reddening, sta r-formation history and 
stellar mass (Sob ral et aLll2010D . The stellar masses of 
these galaxies range fr om ifj 9 - 7-11 - M Q (Table 1; see also 
ISwinbank et al.ll2012D . 

Since the stellar masses are calculated from 2" aper- 
ture photometry (and then c orrected to total m agnitudes 
using aperture corrections, So bral et al.ll2010T) . to esti- 
mate the stellar surface density in the same area as our 
IFU observations, we assume that stellar light follows 
an exponential profile with Sersic index, n scrc = 1 - 2 and 
calculate the fraction of the total stellar mass within 
the disk radius, R (which we define as two times the 
Ha half light radius, r^). Allowing a range of power- 



For local, star-forming galaxies, the exponent, 
n ~ 1.5 and the a bsolute sta r forma tion efficiency, 
v4 = 1.5±0.4x 10~ 4 (|Kennicuttl H998aft implying an 
efficiency for star formation per unit mass of ~ 0.04 
which holds across at least four orders of magnitude in 
gas surface density. 

Combining these relations, the velocity dispersion, er, 
should therefore scale as 



kms 




where Esfr arl d E* are measured in M Q yr 
Mq kpc -2 respectively, R is in kpc, v max in kms" 
G = 4.302 x lO^kpcMg 1 (kms" 1 ) 2 . With a power law 
index of n = 1.4, and a marginally stable disk (Q = l), 
for each galaxy we the refore expect a power law relation 
]gp R + constant (|I 



a cc E°p R + constant (Krumholz et al. 2012). 
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Table 1: Targets fc Galaxy Properties 

ID RA Dec z Ha SFR- Q r* /a <x& Q t& ym E (B - V) log(^) log(^) 

(J2000) (J2000) (M Q /yr) (kpc) (kms" 1 ) (kms" 1 ) 



SHiZELS-1 


0218 26.3 


-04 4701.6 


0.8425 


2 


1.8 + 0.3 


98+15 


112+11 


0.4 + 0.1 


10.03 + 0.15 


9.4 + 0.4 


SHiZELS-4 


10 0155.3 


+02 14 02.6 


0.8317 


1 


1.4 + 0.5 


77+20 




0.0 + 0.2 


9.74 + 0.12 


8.9 + 0.4 


SHiZELS-7 


0217 00.4 


-05 0150.8 


1.4550 


8 


3.7 + 0.2 


75+11 


145+10 


0.2 + 0.2 


9.81 + 0.28 


9.8 + 0.4 


SHiZELS-8 


0218 21.0 


-0519 07.8 


1.4608 


7 


3.1 + 0.3 


69+10 


160+12 


0.2 + 0.2 


10.32 + 0.28 


9.8 + 0.4 


SHiZELS-9 


021713.0 


-04 54 40.7 


1.4625 


6 


4.1 + 0.2 


62+11 


190+20 


0.2 + 0.2 


10.08 + 0.28 


9.8 + 0.4 


SHiZELS-10 


0217 39.0 


-04 44 43.1 


1.4471 


10 


2.3 + 0.2 


64+8 


30+12 


0.3 + 0.2 


9.42 + 0.33 


9.9 + 0.4 


SHiZELS-11 


021821.2 


-05 02 48.9 


1.4858 


8 


1.3 + 0.4 


190+18 


224+15 


0.5 + 0.2 


11.01 + 0.24 


10.1 + 0.4 


SHiZELS-12 


021901.4 


-04 5814.6 


1.4676 


5 


0.9 + 0.5 


115+10 




0.3 + 0.2 


10.59 + 0.30 


9.6 + 0.4 


SHiZELS-14 


10 00 51.6 


+02:33 34.5 


2.2418 


27 


4.6 + 0.4 


131+17 




0.4 + 0.1 


10.90 + 0.20 


10.1+0.4 


Median 






1.46 


7 + 2 


2.4 + 0.7 


75+19 


147+31 


0.3 + 0.1 


10.25 + 0.50 


9.8 + 0.2 



Notes: a Ha star formation rate using the calibration from Kennicutt ( 1998a) with a Chabrier IMF; SFRjja =4.6 x 10 42 Ljja ■ b Ha half 
light radius, deconvolved for the PSF. c Average velocity dispersion for each galaxy, corrected for beam-smearing due to the PSF. Dasym 
denotes the best-fit asymptotic rotation speed of the galaxy, and is corrected for inclination (sec Swinbank ct al. 2012 for details on the 
kinematic modeling of these galaxies). 



law index from n = 1.0 -1.8 and a ra tio of stellar- to 
gas- velocity dispersion of / ff = l-2 (|Korchagin et al.l 
2003), we calculate the best-fit absolute star forma- 
tion efficiency, A and in Fig. [2] we overlay the best- 
fit solutions. Over the range n = 1.0 -1.8, the best 
fit absolute star formation efficiency for the sample is 
A= (4.1 ±2.4) xl0~ 4 M Q yr" 1 kpc -2 (where the error- 
bar incorporates the galaxy-to-galaxy variation, a range 
of / CT = 1 2, and the errors on the stellar masses of each 
galaxy). We note that at low star formation rates and 
stellar masses, there is a non-zero velocity dispersion due 
to the sound speed (c s ) of the gas (c s <, lOkms -1 for the 
Milky Way at the solar circle) which we have neglected 
since this is below both the resolution limit of our obser- 
vations and the minimum velocity dispersion caused the 
stellar disks in these systems. 

We can improve these constraints further assuming 
that star formation in each galaxy behaves in a similar 
way. We reiterate that this model assumes the star for- 
mation is occurring in a marginally Toomre stable disk, 
where the star formation follows the Kennicutt-Schmidt 
Law. Over a range A = 10~ 5 - 10~ 2 (M yr _1 kpc -2 ) and 
n = 0.8 -2.5 we construct a likelihood distribution for 
all nine galaxies and then convolve these to provide 
a composite likelihood distribution, and show this in 
Fig. [3] Although the values of n and A are clearly de- 
generate, the best-fit solutions have n = 1.34±0.15 and 
A = 3.4+ 2 ^ x 10" 4 Mq yr -1 kpc~ 2 Our derived values for 
the absolute star-formation efficiency, A, and power-law 
index, n are within the 1-a of the values derived f or local 
galaxies (e.g. lKennic"uttlll998bt iLerov et al.ll2008D . 

Usi ng the 12 CO to trace the cold molecular 
gas, IGenzel et al.l ((2010) showed that gas and 
star-formation surface densities of high-redshift 
(z ~ 1.5) star- forming galaxies and ULIRGs are 
also well described by the Kennicutt-Schmidt 
relation with coefficients n = 1.17±0.10 and 
A = (3.3 ±1.5) xlO" 4 M Q yr~ 1 kpc~ 2 , which is com- 
parable to the coefficients we derive from our sample. 

In Fig.[3]we plot the star formation and gas-surface sur- 
face density for both local and high-redshift st ar-forming 
galaxies and ULIRGs from lGenzel et all (|2010l ) and over- 
lay the range of acceptable solutions implied by our data. 
We reiterate that we have adopted Q = 1 for this analy- 
sis and note that if we adopt Q < 1 then the absolute 
star formation efficiency will be increased proportionally 



(as shown in Fig. [3]). Nevertheless, this shows that the 
values of n and A we derive are consistent with the lo- 
cal and high-redshift star-forming galaxies and ULIRGs, 
but free from uncertainties associated with converting 
12 CO luminosities to molecular gas mass, CO excitation 
or spatial extent of the gas reservoir. 

Using the values of n and A we have derived, 
we infer cold molecular gas masses for the galaxies 
in our sample of M gas = 10 9-10 M Q with a median 
M gas = 7±2x 10 9 M Q . This suggests a cold molecular 
gas fraction of M gas / (M gas +M*) = 0.3±0.1 but with a 
range of 10-75%, similar to those de rived for other high- 
redshift starburst s in other surveys dTacconi et al.ll20ll 
IDaddi et al.ll2010t ISwinbank et al.ll201lD . 

Finally, with estimates of the disk surface density, we 
can use Eq. Q] to construct maps of the spatially resolved 
Toomre parameter, Q(x, y). Since we set Q = 1 to derive 
the coefficients n and A, by construction the average Q 
across the population is unity, but the relative range of 
Q(x, y) within the ISM of each galaxy is unaffected by 
this assumption. In Fig. Q] we show the maps of Q{x,y) 
for each galaxy in our sample (with contours marking 
Q(x, y) = 0.5, 1.0 and 2.0). This shows that there is a 
range of Toomre Q across the ISM, and to highlight the 
variation with radius, in Fig. 0] we show the Toomre pa- 
rameter within each pixel of each galaxy as a function 
of radius (normalised to the half light radius, r^). This 
shows that in the central regions, on average the Toomre 
Q increases by a factor ~4x compared to Q at the half 
light radius, whilst a radii greater than 7"h, Q decreases 
by approximately the same factor. 

3.1.1. Identification of Star- Forming Regions 

As Fig. [1] shows, the galaxies in our sample ex- 
hibit a range of Ha morphologies, from compact 
(e.g. SHiZELSll & 12) to very extended/clumpy (e.g. 
SHiZELS7, 8, 9 & 14). To identify star-forming re- 
gions on ~ kpc scales and measure their basic proper- 
ties we isolate the star-forming clumps above the back- 
ground (<7bg) by first converting the Ha flux map into 
photon counts (accounting for telescope efficiency) and 
then search for 3<7b g over-densities above the radially av- 
eraged background light distribution. In this calculation, 
we demand that any region is at least as large as the 
PSF. We identify eleven such regions and highlight these 
in Fig. [Q 
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Fig. 3. — Left: The likelihood distribution for the power-law index (n) and absolute star formation efficiency (A) in the Kennicutt-Schmidt 
Law derived from the Hgpn-a relations in Fig. l2l and assuming that the galaxies are marginally unstable, Q = l (equation |4)l . The 
best-fit solutions (within the 1 a contour) have n = 1.34 ± 0.15 and A = 3.4j^'| X 10 — 4 Mg yr -1 kpc -2 . The arrows shows how the absolute 
star formation efficiency would change if we adopt Q = 0.5, or Q = 2 (e.g. Lcroy ct al. 2008). We also plot the position of the clumps 
(adopting n= 1.34). Right: The relation between star formation and gas-surface surface density for local- and high-redshift- star- forming 
galaxies and ULIRGs (Genzel ct al. 2010). The dashed line and shaded region shows the Kennicutt-Schmidt relation with our coefficients 
of n = 1.34 ±0.15 and A = 3.4^t 2 'g X 10 -4 Mq yr - 1 kpc -2 . The solid line shows that best-fit solution for the "Universal" relation from 
I Genzel ct al. (2010), which is well matched to our derived values. 



It is still possible that selecting star-forming regions in 
this way may give misleading results due to random asso- 
ciations and signal-to-noise effects. We therefore use the 
Ha surface brightness distribution from the galaxies and 
randomly generate 10 mock images to test how many 
times a "clump" is identified. We find that only 2 ± 1 
spurious clumps (in our sample of eleven galaxies) could 
be random associations. 

Next, we extract the velocity dispersion and luminos- 
ity of each clump from the using an isophote defining 
the star-forming region and report their values in Ta- 
ble 2 (the clump velocity dispersions have been cor- 
rected for the local velocity gradient from the galaxy 
dynamics and sizes are deconvolved for the PSF). Us- 
ing the velocity dispersion and star formation density 
of each clump, and fixing the power-law index in the 
Kennicutt-Schmidt relation to n = 1.34, we compute 
their absolute star formation efficiencies, deriving a me- 
dian A c i ump = 5.4 ± 1.5 x 10~ 4 (Fig. El). This corresponds 
to an offset (at fixed n) from the galaxy-average of 
^4-ciump / A — 1.3 ± 0.4. Equivalently, if we fix the ab- 
solute star formation efficiency to that of the galaxy- 
average, then the Toomre parameter in these regions is 
Q = 0.8±0.4. 

3.2. The Scaling Relations of Local and High-Redshift 
Star-Forming Regions 

The internal kinematics and luminosities of Hn re- 
gions in local galaxies, derived from the line widths 
of their emission lines, have been the subject of var- 
ious s tudies for some time (e.g. ITerlevich fc Melnickl 
19811: lArsenault et all [19901 iRozas et all 119981 1200a 



Table 2: Physical Properties of the Star-Forming Clumps 



Galaxy 


SFR 


"Ho 


[Nil] /Ha 


rh 




(Moyr- 1 ) 


(km s —1 ) 




(kpc) 


SHiZELS-7 


0.5±0.1 


40 ±10 


0.07 ±0.03 


0.8 ±0.2 


SHiZELS-7 


1.3 ±0.1 


61 ±12 


0.34 ±0.03 


1.0±0.2 


SHiZELS-8 


2.0 ±0.1 


79 ±10 


0.36 ±0.03 


0.7 ±0.2 


SHiZELS-8 


1.6±0.2 


95 ±14 


0.26 ±0.04 


0.8 ±0.2 


SHiZELS-8 


1.9 ±0.1 


140 ± 20 


0.21 ±0.04 


0.9 ±0.2 


SHiZELS-9 


2.1±0.2 


97 ±15 


0.31 ±0.04 


0.7 ±0.2 


SHiZELS-9 


2.3±0.1 


80 ±10 


0.26 ±0.03 


1.3±0.2 


SHiZELS-9 


0.9±0.1 


86 ±14 


0.40 ±0.03 


<0.7 


SHiZELS-14 


0.5±0.1 


56 ±12 


0.12 ±0.04 


0.9 ±0.2 


SHiZELS-14 


1.1±0.2 


121 ±20 


0.24 ±0.03 


<0.7 


SHiZELS-14 


0.2±0.1 


100 ± 25 


-0.03 ±0.05 


0.9 ±0.3 


Median 


1.4 ±0.4 


88 ±9 


0.24 ±0.06 


0.85 ±0.10 



Relano et al.ll2005l ). In particular, if the large line widths 



of star-forming Hii regions reflect the virialization of the 
gas then they can be used to determine their masses. 
However, it is unlikely that this condition holds exactly 
at any time during the evolution of a Hn region due 
to the input of radiative and mecha nical energy, princ i- 
pally from their ionizing stars (e.g. iCastor et al.lll975l) . 



Notes: Half light radius, r^, is deconvolved for PSF and 
the ve locit y dispersion, <r, is corrected for local velocity gradient 
(see § 13.111 . The star formation rates (SFR) are calculated from 
the Ha line luminosity using SFR Hq = 4.6 x 10" 42 Ljja- 

Nonetheless, the least evolved Hn regions may well be 
within a factor of a few (2-3) of having their kinematics 
determined by their virial masses (at an early stage, the 
stellar ionizing luminosities are maximiz ed whereas the 
mech anical energy input is minimized; iLeitherer et al.1 
fl999h . In the case of Hn regions close to virial equi- 
librium, the use of the line-width to compute gaseous 
masses offers a relatively direct means to study the prop- 
erties since it is independent of the small-scale structure 
( density, filling factor etc.). 

ITerlevich fc Melnickl (|1981[ ) showed that the H/3 lumi- 
nosity of the most luminous Hn regions varies as L(H/3) oc 
cr 4 0±0 8 . This result suggests that the most luminous Hn 
regions are likely to be virialized, so that information 
about their masses, and the resultant mass-luminosity 
relation, could be obtained using the virial theorem (they 
also claimed a relation between a radius parameter and 
the square of the velocity dispersion a for Hn regions, 
as further evidence for virial ization). However, more re- 
cent studies, in particular by IRozas et al.l (|2006|) suggest 
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that in super-giant Hn regions, Lex cr 2 - 9±0 - 2 may be a 
more appropriate scaling (the lower exponent arises since 
Hll regions with the largest luminosities are generally 
density-bound, which means that a significant fraction 
of the ionizing radiation escapes and so does not con- 
tribute to the luminosity, making shallower slopes phys- 
ically possible). 

To investigate the scaling relations of star-forming 
regions, in Fig. [5] we show the relations between lu- 
minosity, size and velocity dispersion of the clumps 
in our sample compared to Giant Molecular Clouds 
(GMCs) and Hn regions in the Milky Way and lo- 
cal galaxies (iTerlevich fe Melnickl 119811: Arsenault et alj 



cal galaxies ( lerlcvicn &j lvielmck iy»l; Arsenault et al. 
1990: iBordalo fe TeUed 120111: iFuentes-Masip et all 120001: 



Rozas et al. I 120061) . In this plot, we also include 



the measurements of giant star-forming regions from 
other high-redsh i ft sta r-forming galaxies at z ~ 1 from 
IWisnioski et al.l (I2012D. t he z ~ 1 - 2 galaxies from 
SINS (|Genzel et all ' T2011D . and the clum ps identified 
in str ongly lensed z ~ 1.5-3 galaxies from Uones et"aH 
pOloh and lStark etafl ((20081) . 

Despite the scatter, the radius — a and a - Luminosity 
relations of the high-redshift clumps approximately fol- 
low the same scaling relations as those locally, but ex- 
tending up to ~ kpc scales. Indeed, including all of the 
data-points in the fits, we derive the scaling between size 
(r), luminosity (L) and velocity dispersion (a) of 

log (J_)= (1.01 ±0.08) log (j^r)+(0.8 ±0.1) 

(5) 



and 



log 



L 



ergs" 



= (3.81 ±0.29) log ( — — r ) +(34.7±0.4) 
V km s 1 / 



(6) 

Equation [5] suggests a cx R. If the clouds are self- 
gravitating clouds with a oc i?, then the virial density 
is constant. The relation Loc g 3 81= l=cl 29 is in reasonable 
agreem ent with the early work from lTerlevich fc Melnickl 
(1981), and steeper tha n that found for s uper-giant Hn 
regions in local galaxies (|Rozas et al.ll2006f ) , although the 
large error bars (on both the local and high-redshift data) 
preclude any firm conclusions. Clearly a larger sample is 
required to confirm this result and/or test whether the 
scatter in the data is intrinsic. 

If the star-forming regions we have identified are short 
lived, then these scaling relations effectively reflect ini- 
tial collapse conditions of the clump as it formed, since 
a clu mp can not evolve fa r from those initial conditions 
(e.g. iCeverino et al.|[2010() . In this case, the relation be- 
tween radius, velocity dispers ion and gas mass should fol- 
low r = er 2 / (7rGSdi S k) (see § I3.3|) . In Fig. [5] we therefore 
overlay contours of constant gas mass in the r— a plane, 
which suggests that the initial gas masses for the clumps 
is M™ tM/ = 2 ± 1 x 10 9 M Q a factor - 1000 x more mas- 
sive then the star-forming complexes in local galaxies 
(e.g. see also lElmegreen et al.l l2009t IGenzel et al.l 1201 ll : 
IWisni oski et al.l l2012t ). Assuming our gas mass estimates 
from § 13. 1| then these star-forming regions may contain 
as much as ~10-20% of the cold molecular gas in the 
disk. 

Turning to the relation between size and luminosity of 
the star-forming regions, it is evident from Fig. [5] that 
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Fig. 4. — The variation in Toomre parameter (Q(x, y)) within the 
ISM of the nine galaxies in our sample as a function of (normalised) 
radius. The solid points denote the measurements at each pixel 
within each galaxy and the grey region shows the 18 and 81 %-ile 
limits of the distribution. By construction, the average Toomre Q 
in the sample is Q(x, y) = 1, but varies by a factor ~ 10 within the 
ISM, with the highest-Q (most stable) in the central regions. 

the star formation densities of the high-redshift clumps 
higher than those locally. Indeed, local star-forming re- 
gions follows a scaling relation 



log 



L 



ergs 



(2.91 ±0.15) log 



kpc 



(32.1 ±0.3) 
(7) 

We do not have sufficient number of objects or the dy- 
namic range to measure both the slope and zero-point of 
the size-luminosity relation in the high-redshift clumps, 
and so instead we fix the slope of the local relation (which 
is L oc r 2 - 91 ± 015 ) an d fit for the zero-point evolution and 
obtain 



log 



ergs 



kpc 



(2.91 ± 0.15) log ( — ) + (33.2 ± 0.4) 

This suggests that high-redshift star-forming regions 
have luminosities at a fixed size that are on aver- 
age a factor 15±5x larg er than those locally (see 
also iSwinbank et all [20091 [20lot Uones et al.1 WWt 
IWisnioski et al.ll2012D . We note that high luminosities 
at fixed size have b een found in local s tarbursts, such 
as in the Antennae (jBastian et al J 120 061) . whilst offsets 
of factors ~ 50 x have been inferr ed for star-forming re- 
gions in high-redshift galaxies (e.g. ISwinbank et aP2009t 
I Jones et al.ll2lHl IWisnioski et al.|[2"012l ). 

3.3. The Relation between the Disk and Clump 
Properties 

It is possible to relate the properties of the clumps 
to th e overall properties of the disk (e.g. 1 Ho pkins et al.l 
I2012D . For example, the velocity dispersion of the fastest 
growing Jeans unstable mode which can not be stabilised 
by rotation in a gas disk is given by 



a t {R) 2 = ^GE disk i? 



(9) 



(e.g. lEscala fc Larson! 120081: lElmegreen! 120091: 
iDekel et al.1 120091: IGenzel et all l201lT iLivermore et all 
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Fig. 5. — Scaling relations between size, luminosity and velocity dispersion for the star-forming regions in our high-r edshift galaxies com- 
pared t o those in local GMCs a n d Hll regions. In a l l of th e se plots, we baseline our mea s urements aga i nst loc al data from TerlevicE & Melnickl 
yip); lArsenault et~aTI (fT590T) ; IBordalo fc Telled 1(2011) : IFuentes-Masip et al.l (|2000h : IRozas et ahl j2uM) . Left: The relation between ve - 
locity disp ersion and size. For th e high-redshift star-forming regio ns, we also include clump measurements from SINS (Genzel et al, 120111) , 
ZWigglcs (Wisnioski et al. 2012) and the cluster arc survey from (Jones et al. 2010). We also plot the properties of the HiZ GMCs from 
the numerical simulations from Hopkins (2012). The dashed line shows a fit to the data of the form r oc <r 101 . The dashed lines show 
lines of constant gas mass (|S). Middle: The relation between velocity dispersion and luminosity of star-forming regions in high-redshift 
galaxies compared to those locally. The dashed line denotes Loc cr 3 8 which provides a good match to both the local and high-redshift 
data. Right: The scaling relation between size and luminosity of star-forming regions. The high h igh-redshift star-formi ng regions have 
luminosity densities which are a factor ~15±5 X higher than those typically found locally (see also Wisnioski et al. 2012). 

(appropriate for a uniform density sphere). Although 
this calculation should be considered crude as it is 
unclear whether the velocity dispersions we measure are 
virial, it is encouraging that the predicted surface mass 
densities of the clumps are similar to those inferred from 
their velo city dispersions and size s. 

Fi nally. IHopkinsI (120121) ( see also lEscala fc Larsonll2008l 
and iLivermore et all |2012[ ) show that for a marginally 
stable disk of finite thickness, density structures on scales 
greater than h will tend to be stabilised by rotation which 
leads to an exponential cut off of the clump mass function 
above 



I2012D . The critical density for collapse (p c ), on scale R 
from a turbulent ISM is given by 



Pc 



4vri? 3 



M, 



8irR 2 G 



<Jt{R? 



(10) 



where at (R) is the line of sight velocity turbulent velocity 
dispersion and Mj is the Jeans mass. The critical density 
for collapse therefore scales as 



Pc(R) 



9 y 



(11) 



Assuming that the cloud contracts by a factor ~ 2.5 as it 
collapses, the post-collapse surface density of the cloud 
is 

£ c , oud ~ 10 Pc R~ 10E disk (12) 

(see also ILivermore et"al~ll2012[ ). Thus, the surface den- 
sity of the collapsed cloud is independent of radius and 
proportional to the surface density of the disk, with the 
normalisation set by the collapse factor and under the 
assumption Q = l. IHopkinsI (|2012f) show that this model 
provides an reasonable fit to giant molecular clouds in 
the Milky- Way, and further, suggests that the surface 
density (and hence surface brightness) of clouds should 
increase with the surface density of the disk. 

Using our estimates of the stellar and gas 
masses and spatial extent of the galaxies in 
our sample, we derive disk surface densities of 
Sdisk = 1-1 ± 0.4 x 10 9 M Q kpc -2 , and hence expect 
the mass surface densities of the star-forming re- 
gions that form to have mass surface densities of 
^ciump ~ 10 10 M© kpc -2 . It is instructive to compare 
this to the average mass surface density of the clumps. 
For example, assuming that their velocity dispersions 
are virial and adopting M c \ ump = Ca 2 rh / G, using 
the average velocity dispersion and size of the clumps 
(Table 2), we derive an average clump mass surface 
density of S c ium P = 8 ± 2 x 10 9 M Q kpc" 2 with C = 5 



f Pc(h)h* 



3^G 2 EL k 



(13) 



M 

M Q 



= 8.6 x W 



isk 



10 M pc" 



100 kms -1 kpc 



(14) 

This suggests that the most massive clumps that can 
form in a disk ( "the cut off mass" ) depends strongly on 
the disk surface density - increasing the disk surface den- 
sity i ncreases mas s of the clumps that are able to form 
(e.g. [Escala 2 01 lh . However, there is also a competing 
(stabilising) factor from the epicyclic frequency such that 
a fixed radius, higher circular velocities reduce the mass 
of the largest clumps able to form. 

Applying equation [14] to the Milky- Way, with a cold 
mole cular gas fraction of 10%, /cr = 2 (|Korchagin et al.l 
2003), the average surface density is Edisk = 35 M Q pc 
and for n = 220 km s" 1 / 8 kpc (|Feast fc Whitelockl[T997t l 
the cut off mass should be Mo ~ 10 7 M Q , in good agree- 
ment with t he characteristic m ass of the largest galactic 
GMCs (e.g. IStark fc Ledl200fil ) 

How does the cut-off mass for our high-redshift sam- 
ple compare to local galaxies? For f a = 2, and us- 
ing the scaling relations derived in § 13.11 to estimate 
the gas mass (Table 1), (A = 3.4 x 10" 4 M yr" 1 kpc -2 
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and n = 1.34) we derive a range of cut off masses of 
Mo = 0.3-30 x 10 9 M Q (with a median and error of the 
sample of M = 9 ± 5 x 10 9 M Q ). This is similar to 
the mass inferred for the brightest star-forming regions 
seen in high-resolution images of other high-r e dshift 
galaxies (lElmegreenlll989t lElmegreen et all 120071 I200S ' 



Bournaud & Elmcgreen 2009; IForster Schrcibcr et al 



2011bt iGenzel et all 120111: iWisnioski et all I2012D . and 



a factor ~ 1000 x higher than the largest characteristic 
mass of a star-forming region in the Milky- Way. 

In Fig. |5] we plot our estimates of the the cut off mass 
versus the clump star- form ation densities for the galax- 
ies in our sample (see also iLivermore et all [20121 We 
use the Ha derived star-formation rate for each clump, 
corrected for galaxy reddening (note that we do not have 
reddening estimates for individual clumps and so we as- 
sume a factor 2 x uncertainty in their star formation sur- 
face density). We also include estimates of the cut off 
mass and star formation surf ace density fr o m the SINS 
survey of z ~ 2 galaxies fr om IGenzel et all (120111) (with 
dyn amics measured fr om IForster Schreiber et all 120091 
and iCresci et all 120091 ) , as well as mea s ureme nts from 
the lensing samples o f ILivermore et all (|2012l ) [z ~ 1) 
and lJones et~a l. (2010) (z~2). Although the error-bars 
on individual measurements are large (particularly due 
to the uncertainties in deriving the gas surface density 
from the Kennicutt-Schmidt law), as can be seen from 
Fig. [51 galaxies with high cut-off masses tend to have 
higher clump luminosity surface densities. 

It is also useful to adopt simple models for the evolu- 
tion of galaxy disks and gas fraction to investigate how 
the cut off mass and clump properti es may be expe c ted to 
evolve with redshift. For example, iDutton et all ([201 ID 
present an analytic model for the evolution of disk scal- 
ing relations (siz e, rotational velocit y and stellar mass 
with redshift; sec IDutton et ai][20TTI Table 3). Combin- 
ing with a simple model for the evolu tion of the gas frac- 
tion / gas cx(l + z) b with 6 = 1.5-2.5 ([Geach et al.ll201lD 
and using Eq. [12] and Eq. [14] we show the expected evo- 
lution of the cut-off mass and clump luminosity surface 
density with redshift. This shows that as the gas fraction 
increases (and adopting evolving models for the size, disk 
and circular velocity of galaxies), then the cut off mass 
should increase by a factor 10- 100, x over the redshift 
range z = 0-2.5 whilst the star formation density of the 
clumps should increase by a approximately an order of 
magnitude over the same redshift range. Although this is 
a simple model, this framework allows us to understand 
why the properties of the star-forming clumps within the 
ISM of our sample of high-redshift galaxies are different 
to those typically found in star-forming galaxies locally. 

4. CONCLUSIONS 

We have presented resolved spectroscopy of nine star- 
forming galaxies at z = 0.84-2.23 selected from the 
UKIRT/HiZELS survey. These galaxies have reddenning 
corrected star-formation rates of SFR= 16±5M Q yr _1 
and so are representa tive of the high-redshift population 
(jSobral et all l2012bh . The Ha dynamics suggest that 
the ionised gas in at least six galaxies is in the form of 
large, rotating disks. We use the inferred rotation speeds 
of these systems, together with the spatial extent of the 
Ha to investigate the star formation within the ISM, and 
we derive the following main conclusions: 
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Fig. 6. — The most massive clumps that can form (the "cut off 
mass", Mo) as a function of clump star formation surface den- 
sity for SHiZELS galaxies. The cut off mass is related to the 
disk surface density (Edisk) and epicyclic frequency (k) via Mo oc 
S^ isk K — 4 . The z = observations are derived from the Spitzer In- 
frared Nearby Galaxy Survey ( SINGS) Hvennicutt et all J2T)03T) . We 
also include in the plot measurements of other high-re dshift star- 
forming galaxies from the SINS s u rvey jG cnzcl et al. 201 J) and 
lensing surveys (Jones et al. 2010; Livcrmorc et al. 2012). This 
shows that the cut off mass and star-formation surface densi- 
ties of the high-redshift star-forming regions are (up to) a fac- 
tor ~ 100 X higher than star-forming regions in local galaxies. 
Using a simple model for galaxies with evolving gas fractions 
(/gas <x (1 + z)' 2 ^ °' 5 ' ) and using the redshift evolution of disk 
scaling relations (si ze, rotational velocity an d ste llar mass) from 
IDutton et al. (2011.) and using equation 1121 and 1141 can be used 
to derive model tracks to show how the cut off mass and clump 
star formation surface density are expected to evolve with redshift. 
We plot these tracks for a ratio of stellar-to-gas velocity dispersion, 
fcr = 2, but also show how the results change if we instead adopted 
vary for fa- = 1. These tracks shows that the cut off mass and 
clump star formation surface density should increase by 1 — 2dex 
between z = and z = 3. 

• The star formation and velocity dispersion within 
the ISM of these high-redshift galaxies follow a power- 
law relation of the form a oc ASga" + constant where the 
coefficients, A and n are set by the Kennicutt-Schmidt 
Law (Ssfr = J^Ilgas) and the constant includes the disk 
stellar surface density of finite thickness. Assuming the 
gas disks are marginally stable (Q = 1) we combine the 
solutions for each galaxy and derive best-fit parameters 
of power-law exponent, n= 1.34±0.15 and absolute star 
formation efficiency, A — 3.4^ 2 g x 10~ 4 M Q yr" 1 kpc~ 2 . 
These values are consistent with the parameters derived 
via 12 CO observations for both local and high-redshift 
star-forming galaxies, but free from any assumptions 
about 12 CO-H2 conversion factors, 12 CO excitation or 
the spatial extent of the gas reservoir. 

• Applying these coefficients, we infer cold molecular 



gas - 



-TO 9 - 10 



Mq with a me- 



gas masses in the range M 

dian M gas = 7 ± 2 x 10 9 M Q and hence a cold molecular 
gas fraction of M gas / (M gas + M*) = 0.3 ± 0.1 but with a 
range of 10-75%. 

• Using a simple analytic model, we show that the 
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largest structures that can form within the disk (the cut- 
off mass, Mo) are set by the disk surface density with a 
competing (stabilising) force from the epicyclic frequency 
such that Mo oc £^ isk K . For the galaxies in our sam- 
ple, we derive cut off masses of Mq ~ 10 9 M©, a factor 
~ 1000 x higher than the largest characteristic mass of 
GMCs in the Milky- Way. 

• Within the ISM of these galaxies, we reliably iso- 
late eleven ~ kpc-scale star-forming regions and measure 
their properites. We show that their luminosities and 
velocity dispersions follow the same scaling relations be- 
tween size and velocity dispersion as local Hn regions. 
Assuming the line widths are virial, the masses derived 
for these star-forming regions are consistent with those 
implied by the cut-off mass. However, we find that the 
luminosity densities of these star-forming regions are a 
factor ^15x higher than those typically found locally, 
which we attribute to the requirement that the surface 
density of the (collapsed) cloud must be ~ 10 x that of 
the disk. 

Overall, the scaling relations we have derived suggest 
that the star formation processes in high-redshift disks 
are similar to those in local spiral galaxies, but occurring 
in systems with a gas rich and turbulent ISM. Given the 
pau city of gas-rich, clumpy disk-like high-redshift galax- 
ies (lElmegreen et al.ll2007ll2009tlForster Schreiber et al.l 
1201 lal) . the next step in these studies is to spatially re- 
solve the cold molecular gas via CO spectroscopy in a well 
selected sample in order to better constrain the intera- 
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tion between star-formation and gas dynamics. Through 
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lations (e.g. iCrain et al.ll2009t Ivan de Voort et al.ll2011t 
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to differentiate whether the dominant mode of accretion 
is via three-dimensional cold gas flows accrete from the 
inter-galactic medium, or from two-dimensions from out- 
skirts of the disk as gas cools from the hot halo. 
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